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ABSTRACT 

Using a dynamical 3-D reconstruction procedure we estimate the peculiar ve- 
locities of i? > Abell/ACO galaxy clusters from their measured redshift within 
25000 km/sec. The reconstruction algorithm relies on the linear gravitational insta- 
bility hypothesis, assumes linear biasing and requires an input value of the cluster 
/3-parameter (/3c = 0°'^/6c), which we estimated in Branchini & Plionis (1995) to 
be l3c — 0.21. The resulting cluster velocity field is dominated by a large scale 
streaming motion along the Perseus Pisces-Great Attractor base-line directed to- 
wards the Shapley concentration, in qualitative agreement with the galaxy velocity 
field on smaller scales. Fitting the predicted cluster peculiar velocities to a dipole 
term, in the local group frame and within a distance of ~ 18000 km/sec, we recover 
extremely well both the local group velocity and direction, in disagreement with the 
Lauer & Postman (1994) observation. However, we find a ~ 6% probability that 
their observed velocity field could be a realization of our corresponding one, if the 
latter is convolved with their large distance dependent errors. Our predicted cluster 
bulk velocity amplitude agrees well with that deduced by the POTENT and the da 
Costa et al. (1995) analyses of observed galaxy motions at ~ 5000 — 6000 km/sec; it 
decreases thereafter while at the Lauer & Postman limiting depth (~ 15000 km/sec) 
its amplitude is ~ 150 km/sec, in comfortable agreement with most cosmological 
models. 

Subject headings: cosmology: theory - galaxies: distances and redshifts - large-scale 
structure of universe - galaxies: clusters: general 
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1. Introduction 

There is strong observational evidence for the 
existence of coherent large-scale galaxy flows in the 
local universe, extending from the Perseus-Pisces 
region on the one side to the Hydra-Centaurus/Great 
Attractor region on the other, pointing within 
^ 40° of the CMB dipole direction (see reviews 
and references in Dekel 1994 and Strauss & Willick 
1995). These results, together with the large-scale 
coherence of the cluster gravitational acceleration 
(indicated by the fact that the differential clus- 
ter dipole in large equal-volume shells is aligned 
with the CMB dipole in each shell out to ~ 16000 
km/sec; Plionis & Valdarnini 1991; Plionis, Coles 
& Catelan 1993; Plionis 1995), present a consistent 
picture in which the local group [hereafter LG] par- 
ticipates in a large-scale bulk motion induced by 
gravity, encompassing a volume of radius ^ 15000 
km/sec. 

This picture has recently been challenged by 
Lauer & Postman (1994) [hereafter LP94] who 
have extended the cosmic flow studies to very large 
scales using the brightest cluster galaxies as stan- 
dard candles. They find that the LG motion with 
respect to the frame defined by the Abell/ACO 
clusters within 15000 km/sec moves in a direction 
~ 80° away from that of the CMB dipole, which 
then implies that, if the CMB dipole is a Doppler 
effect, the whole cluster frame is moving with re- 
spect to the CMB rest-frame with ~ 700 km/sec. 
Such a large velocity on such large scales is difficult 
to reconcile with galaxy bulk velocities on smaller 
scales, and with the current models of structure 
formation (cf. Strauss et al. 1995; Feldman & 
Watkins 1994). Furthermore it is difficult to un- 
derstand why the LG peculiar acceleration, esti- 
mated from the cluster dipole, is well aligned with 
the mass (CMB) dipole (see Branchini & Plionis 



1995, and references therein), while the observed 
LG peculiar velocity (as estimated by LP94), with 
respect to the same sample of clusters, is not. 

The main aim of this Letter is to predict the 
Abell/ACO cluster velocity field and bulk fiow, 
within the gravitational instability and linear bias- 
ing framework and compare it with determinations 
based on the observed Tpeculiai velocities of galaxies 
(cf. POTENT) and of the LP94 clusters (see also 
Scaramella 1995 for complementary approach). 

2. Method 

We use the linear gravitational instability [GI 
hereafter] framework, linear biasing and an itera- 
tive technique, similar to that of Yahil et al. (1991), 
to reconstruct the Abell/ACO real cluster positions 
within 25000 km/sec, starting from their redshift 
space distribution, and thus obtain their peculiar 
velocity field. In the linear approximation the pe- 
culiar velocity u at position r is proportional to the 
gravitational acceleration (cf. Peebles 1980): 

= tj ^cWpf^d^r (1) 

where /3c = il'^ '^/hc, SIq is the present value of the 
cosmological density parameter and he is the usual 
linear bias parameter that relates the cluster and 
mass overdensities; 6c = h^bm- Note that the large 
relative separations of galaxy clusters causes them 
to sparsely trace the underlying density field; it 
is however this fact that make them ideal probes 
of the linear cosmic dynamics (cf. Bahcall et al. 
1994), although non-linear effects could be present 
in high density regions, ie., up to scales of ~ 1000 
km/sec (cf. Croft k Efstathiou 1994). 

Our technique is a two step procedure whose 
complete description can be found in Branchini & 
Plionis (1995) [BP95 hereafter]; here we just recall 
the general idea. 
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• The first step consists in reconstructing the 
whole sky redshift-space cluster distribution 
with cz < 25000 km/sec. This is done 
by Montecarlo generating a population of 
synthetic objects outside the zone of avoid- 
ance (i.e. at l&l > 20°) and within 20000 
km/s, whose spatial distribution accounts for 
galactic absorption, the radial selection func- 
tion and density inhomogeneities between the 
Abell and AGO catalogues. Furthermore, the 
synthetic clusters are spatially correlated to 
the real ones according to the observed Ccc(f). 
The object distribution within the zone of 
avoidance [ZoA hereafter] is then recovered 
by randomly cloning the cluster distribution 
within redshift-galactic longitude bins in the 
nearby latitude strips. The synthetic clus- 
ters within the [20000, 25000] km/sec shell 
are also distributed in order to minimize the 
known observational biases but spatially un- 
correlated with real clusters, the latter to 
avoid introducing shot noise effects related 
to the exponential decrease of the real cluster 
number density. Finally, the density field be- 
yond 25000 km/s is considered homogeneous 
and isotropic. 

• The second step consists in applying to the 
whole-sky distribution of real -|- synthetic 
clusters reconstruction procedure that itera- 
tively minimizes the reshift space distortions 
allowing to recover the 3-D cluster positions 
and peculiar velocities. This procedure as- 
sumes linear GI theory, linear biasing and re- 
quires an input /3c parameter, which we take 
it to be l3c = 0.21(±0.03), obtained by com- 
paring the amplitudes of the GMB and 3-D 
cluster dipoles (see BP95). To deal with pos- 
sible non-linear effects on small scales we ap- 
ply a top-hat smoothing of the forces. 



An estimate of the intrinsic and possible system- 
atic errors, introduced by modelling the observa- 
tional biases (galactic absorption, radial selection 
and Abell/ AGO homogenization scheme) which we 
call observational error, can be obtained by vary- 
ing the model parameters in a plausible range and 
by generating different realizations of the cluster 
velocity field (for more details see BP95 and Plio- 
nis et al. m preparation). It turns out that for 
the input /3c parameter used and for a smooth- 
ing radius of 10^ km/sec, the observational error 
distribution has a mean and standard deviation 
of ((To) 80(±60) km/sec. Similarly, the intrin- 
sic uncertainty in the reconstruction procedure, (T/, 
obtained as the scatter of the cluster velocities re- 
sulting from different Montecarlo realizations of the 
same model, is (cr/) 144(±78) km/sec. 

Another source of uncertainty is the approxi- 
mate nature of the ZoA model and the possible sys- 
tematic effects related to its increase with distance. 
In BP95 we implemented various schemes for fill- 
ing the ZoA and found that the error on the re- 
constructed cluster positions was smaller than the 
intrinsic one and, more importantly, the resulting 
cluster dipole was nearly unaffected, which implies 
the stability of the bulk velocity measurement. 

The parameters that affect mostly the resulting 
velocity field are the value of the input /3c param- 
eter which is required to start the iterative proce- 
dure and to a lesser extent the smoothing radius, 
Rsm, which affects however the velocities only in 
the high density regions. In BP95 we have verified 
that the choice of the input /3c parameter does not 
bias the reconstruction procedure and that the so 
called "Kaiser effect" is unimportant in the present 
analysis. 
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3. The Local Group Velocity with respect 
to the Clusters 

Similarly to LP94 we solve for the LG peculiar 
velocity with respect to the cluster frame, ulg, by 
minimizing: 

^.^^flEf^^id^u^l^^ (2) 

where di is the reconstructed 3-D cluster distance 
in the LG frame, r,- the unit position vector and 
(Tt j- is the individual cluster total error computed 
by adding in quadrature crji, (Tgi and (T^(= 300 
km/sec) which represents the average uncertainty 
in czi. Note, that the significance measure of 
the dipole fit is ill-defined since the velocity errors 
are coupled. We have, however, investigated the 
stability of our solution to variations of the sample 
size (with > 50% reduction) and error-weighting 
and found that indeed it is very robust. 

In table 1 we present the solution of eq.(2) for 
the Rsm = 1 and 2x10^ km/sec cases and for 
three limiting radii of the volume used (boldface 
quantities refer to our standard case). It is im- 
portant to note that out to the limiting depth of 
the LP94 sample (r^p ~ 15000 km/sec) we find 
|u|lg(''lp) ~ 510 ± 100 km/sec with a misalign- 
ment angle with respect to the CMB dipole apex 
of only 69cmb ~ 5° while its asymptotic value 
(Ri 635 ±70 km/sec) is reached at ~ 18000 km/sec. 
This result, which assumes the GI hypothesis and 
agrees well with the observed CMB dipole, dis- 
agrees with the observed LG peculiar velocity, as 
measured by LP94, which has a similar amplitude 
but 69cmb ~ 80°. This apparent discrepancy could 
be, among other things, due to the large distance 
dependent uncertainties of the LP94 velocities; in 
other words their velocity field could be a realiza- 
tion of an underlying field, represented by our re- 
construction once convolved with the large distance 



dependent errors {a 0.16 r). 

Using Montecarlo simulations, in which we re- 
placed each of our cluster 3-D distances, ri, with 
(1 -|- ii)ri (where each ii is drawn from a Gaussian 
with jj, = Q and a = 0.16) and then fitting vllg [via 
eq.(3)] for the resulting velocity fields, while using 
a z~'^ weighting (see LP94), we have found that in 
~ 6% of the cases the derived vllg was within Icr, 
in amplitude and direction, of that of LP94 (more 
details in Plionis et al m preparation). 

Figure 1 shows |u|lg, the solution of eq.(2), 
computed within spheres of increasing radius to- 
gether with the analogous quantity, |u|g, obtained 
from the gravitational acceleration induced by clus- 
ters (cluster dipole; see BP95) acting on the LG, 
after including a contribution of a ~ 170 km/sec 
Virgocentric infall and for [3^ = 0.21 (error bars 
indicate 1 <tt)- There is a very good matching be- 
tween the two quantities while their misalignment 
with the CMB dipole direction is ^ 12° at the con- 
vergence depth. Although vllg and u^ are not in- 
dependent measures of the dipole, since they de- 
rive from the same underlying density field, their 
good matching constitutes a non-trivial demand- 
ing test of our reconstruction procedure, which we 
pass with success. 

4. The Cluster Velocity Field and Bulk 
Flow 

In figure 2 we present, for the Rsm = 10^ km/sec 
case, the cluster peculiar velocity field in a 8000 
km/sec wide slice projected onto the supergalac- 
tic plane, where most prominent superclusters lie 
(Hydra-Centaurus, Coma, Shapley, Perseus-Pisces, 
Ursa-Major and Grus-Indus). Open and filled dots 
refer to infiowing and outfiowing objects, respec- 
tively, while the length of each vector is equal to 
3 times the line of sight component of the pecu- 
liar velocity in the CMB frame. The small circle 
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at the center represents the typical region spanned 
by dynamical analyses based on galaxy peculiar 
velocities (cf. Dekel 1994 and Strauss & Willick 
1995). The most prominent feature is a large 
coherent motion in the general direction of the 
CMB dipole towards the Shapley Concentration 
{X,Y) = (-13000, +9000) km/sec which does not 
have, however, a constant amplitude; it is small in 
the Perseus-Pisces region (X,Y)= (+8000, -4000) 
km/sec, then rises in the Great Attractor region 
(X,Y)= (—4000, +500) km/sec, while dropping on 
its backside (cf. Dressier & Faber 1991, Mathewson 
et al. 1994). Moreover, the bulk velocity rises again 
near the Shapley concentration where a strong back 
infall is apparent. Other features include the neg- 
ligible peculiar velocities in the Coma region (cf. 
Courteau 1992) and the local infall in the Ursa- 
Major region, at (X,Y)=(7000, 14000) km/sec. 

We have measured the Abell/ACO cluster bulk 
velocity which is defined as the center of mass ve- 
locity of a specified region and is given by the in- 
tegral of the cluster peculiar velocities u(x) over 
a selected volume specified by a selection function 
^(x): 

ylo^r^a.) = [ ""^ H^) u(x) dx (3) 
■Jo 

Assuming that clusters trace the mass, we can write 
eq.(3) for our discrete composite cluster sample 
as Vg^'*^ = WiUi/ Wi, where the sums ex- 
tend over both real and synthetic objects and Wi 
is a weight that accounts for cluster masses and 
Abell/ACO sample differences (however, similar 
results are obtained even for Wi = 1). The dis- 
crete approximation need not to be equal to that 
of eq.(3) since it is a sum over the observed clus- 
ters which are inhomogeneously distributed. In our 
case, however, this sum extends over the whole 
real+synthetic cluster distribution and not only 
over the real objects. Moreover, we have found con- 



sistent estimates of Vg^'*^ using the cluster peculiar 
velocities either at the positions of the clusters or 
at a regular grid with grid size of 2000 km/sec. 

The cluster bulk velocity can be also defined as 
the residual velocity of the whole cluster frame. 
Therefore an alternative estimator of the bulk ve- 
locity, which utilizes only real clusters and only the 
line of sight component of their peculiar velocities, 
which was also used by LP94, is: 

yloHrrna.) = C - ULGirma.) (4) 

where C is the CMB dipole vector and iiLG(fmax) 
is given by eq.(2). Note that our reconstruction 
procedure assumes that the mass fiuctuations re- 
sponsible for the cluster peculiar motions are con- 
tained within the sample considered, imposing the 
bulk fiow to vanish beyond 25000 km/sec. We have 
verified, varying the limiting sample depth and us- 
ing an alternative reconstruction scheme that in- 
creases the amount of clustering in the 20000 - 
25000 km/sec range (see BP95), that this con- 
straint does not appreciably affect the bulk fiow 
amplitude within 20000 km/sec. 

In table 2 we present the results of eq.(4) as a 
function of sample limiting depth for two differ- 
ent smoothing radii. Note that the bulk fiow at 
r^p has an amplitude of ~ 150 km/sec in comfort- 
able agreement with currently accepted cosmolog- 
ical models. In figure 3 we plot, for Rgm = 10^ 
km/sec, both estimators Vg^'*^ and Vj^'*^ as filled 
dots and starred symbols respectively, but now as 
a function of the effective depth (reff ~ 3/4rmaj;). 
Note that for the Vid estimator we have used 
fmax > 10000 km/sec in order to have sufficient 
data for the minimization to be stable. Error 
bars represent 1 a total uncertainties. Although 
the two estimators are not independent, they are 
however based on a different set of clusters and 
velocities and they provide a consistent estimate 
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of the bulk velocity. We also plot as open dots 
the bulk velocity obtained by the POTENT recon- 
struction of the Mark III velocity field (Dekel 1994) 
and as starred symbols the recent da Costa et al. 
(1995) bulk velocity. At the region where both, the 
galaxy and cluster bulk velocity estimates overlap 
(at reff ~ 4000-6000 km/sec ), the difl'erent bulk 
fiow amplitudes appear to be in very good agree- 
ment with each other. Note that our bulk fiow 
determination is mostly unaffected by the window 
'shrinkage' effect (Kaiser 1988), since the cluster 
selection function has a value ~ 1 up to ~20000 
km/sec and does not suffer from the use of only line 
of sight peculiar velocities (Regos & Szalay 1989) 
since the Vg^j*^ estimator is based on the full 3D 
velocity field. 

5. Conclusions 

An iterative reconstruction procedure, based on 
linear GI theory and linear biasing, has been used 
to derive the velocity field traced by Abell/ACO 
clusters within 25000 km/sec. Our main results 
are: 

1. The predicted cluster velocity field, in the LG 
frame, is such that it refiects the whole LG 
motion with respect to the CMB, in apparent 
disagreement with the LP94 velocity data. 
Our derived LG velocity is well aligned with 
the CMB dipole and it reaches its asymptotic 
value at ~ 18000 km/sec in agreement with 
the locally derived cluster dipole. We have 
found, a small but non-negligible probabil- 
ity (~ 6%) that the LP94 observed velocity 
field could be consistent with our predicted 
one; their apparent disagreement being pos- 
sibly due to the convolution of the underlying 
velocity field with their large distance depen- 
dent errors. 



2. The main features of the observed veloc- 
ity field probed by galaxies (cf. Strauss & 
Willick 1995) are reproduced also by our pre- 
dicted cluster velocity field. There is an ev- 
ident extension of our cluster bulk fiow out 
to ~ 15000 km/sec, where a back-infall to 
the Shapley concentration is evident. The 
derived bulk fiow velocity has an amplitude 
in very good agreement with that of PO- 
TENT and da Costa et al. (1995) at ~ 5000 
km/sec, it decreases thereafter while point- 
ing towards the CMB dipole direction. Our 
predicted bulk velocity at rmax ~ 10000 and 
15000 km/sec is ~300 and ~150 km/sec, re- 
spectively, consistent with most theories of 
structure formation. 
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Figure 1: Gomparison between the LG peculiar 
acceleration, |u|g, as derived by the cluster dipole 
and the LG peculiar velocity, |u|lg, derived from 
eq.(2) for the Rgm = 10^ km/sec case. 

Figure 2: The cluster velocity field of a 8000 
km/sec strip projected onto the supergalactic plane 
for the Rsm = 10^ km/sec case. Open dots are 
infiowing clusters while filled dots are outfiowing 
ones. The velocity vectors have been multiplied by 
a factor 3 to enhance their visibility. 

Figure 3: The cluster bulk velocity as estimated 
by the two methods described in the text together 
with the galaxy based POTENT and the da Gosta 
et al. (1995) values. 



7 



Table 1: Local Group motion with respect to the 
cluster frame within distance rmax and for 2 differ- 
ent smoothing radii (both in Mpc). 
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Table 2: Residual bulk motion of cluster frame 
within distance Vmax- 
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